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Abstract 

We evaluate equal time point to point spatial correlation functions of mesonic 
currents at finite temperature. For this purpose we consider the QCD vacuum 
structure in terms of quark antiquark condensates and their fluctuations in 
terms of an irreducible four point structure of the vacuum. The temperature 
dependence of quark condensates is modeled using chiral perturbation theory 
for low temperatures and lattice QCD simulations near the critical temper- 
ature. We first consider the propagation of quarks in a condensate medium 
at finite temperature. We then determine the correlation functions in a hot 
medium. Parameters such as mass, coupling constant and threshold energy 
are deduced from the finite temperature correlators. We find that all of them 
decrease close to the critical temperature. 
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I. INTRODUCTION 



The structure of vacuum in Quantum Chromodynamics (QCD) is one of the most in- 
teresting question in strong interaction physics The evidence for quark and gluon con- 
densates in vacuum is a reflection of its complex nature 0. Determination of correlation 
functions of hadronic currents in such a vacuum state provides rich information re- 
garding interquark interaction as a function of their spatial separation as well as on hadron 
spectroscopy. These are some of the nonperturbative feature of QCD and are of great value 
in understanding the ground state structure of the theory of strong interactions j^Q. 

We have studied mesonic and baryonic current correlators at zero temperature with a 
non-trivial structure for the ground state with quark antiquark condensates 1^,^. It was 
shown that the square of the quark propagator does not reproduce the correlation function 
for the pion deduced from phenomenology. In order to match the data it was necessary 
to introduce an irreducible four point structure for the quarks in the vacuum. This may 
be looked upon as an effective way of incorporating gluon condensate contribution to the 
correlator. 

As is well known the QCD vacuum state changes with temperature. Lattice monte 
carlo simulations suggest that chiral symmetry is restored around 150 MeV. In view of this 
the present note is aimed at looking at the behaviour of the meson correlation functions at 
finite temperature. This is of great interest in the context of behaviour of hadrons around 
the chiral phase transition associated with quark gluon plasma PJ^. It may be noted that 
there is little phenomenological information in this regime but there are several theoretical 
studies essentially using sum rule methods. The main objective here is to employ a 



different nonperturbative approach developed by us ||Tl|. This has been sucessful at zero 
temperature, and, its extension to finite temperature is therefore of interest. In particular, 
we will obtain temperature dependence of masses, coupling constants and threshold energies 
for the pion and rho mesons. 

We organise the paper as follows. In section II we discuss the quark condensate at finite 
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temperature to fix the parameter appearing in the ansatz of the ground state of QCD. We 
then discuss in section III the quark propagation in the thermal vacuum. In section IV we 
calculate meson correlation functions at finite temperature. Finally we discuss the results 
in section V. 



II. QUARK CONDENSATE AT FINITE TEMPERATURE 

To calculate the correlators at finite temperature we need the expression for the equal 
time propagator for the interacting quark field operators. We have developed earlier a 
vacuum structure in terms of quark antiquark condensates with a condensate function h{k). 
The equal time propagator could then be calculated in terms of the condensate function 
0. One can generalise this to finite temperature using the method of thermofield dynamics. 
Here the thermal average is obtained as an expectation value of the operator over the thermal 



vacuum 



T^ . This leads to 



iU^\.lA(n\\^ - J!l- f e-'^-'^A-paik,T)dk (1) 

The thermal vacuum is obtained from the zero temperature vacuum by a thermal Bogoli- 
ubov transformation in an extended Hilbert space involving extra field operators (thermal 
doubling of operators) [0. The functions A± (Eq. |T]) for the case of two flavour massless 
quarks, are given as(with k = \k\), 
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A±(fc,T) 2 



1 ± cos 2^(7° sm2h{k) + a-k cos2h{k)) 



(2) 



In the above h(k) is the condensate function [^^|TT| corresponding to the Bogoliubov trans- 
formation to include a condensate structure in the vacuum. The function 6 is associated 
with the thermal Bogoliubov transformation and is related to the distribution function as 

m 



sin^ e(k) = , (3) 
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(3 being the inverse temperature. Further e{k) is the single particle energy given as e{k) 



kP' + m{ky. In the presence of condensate the dynamical mass is given as m{k) = m + 
kta.n2h{k), m being a possible current quark mass 0. 

We had earlier taken a gaussian ansatz for the condensate function sm2h{k) = e.-^^^^l'^. 
In order to determine the parameter i?, we had taken a value of R consistent with hadronic 
correlator phenomenology. We choose a similar structure for the condensate function at finite 
temperature namely sm2h{k) = ^-^iTf^'^/'^ with R(T) now being temperature dependent. 

In order to determine R{T) or equivalently the ratio S{T) = R{T = 0)/R{T), we first 
evaluate our expression of the order parameter (the condensate value) at finite temperature. 
In terms of the dimensionless variable t] = Rk, this is given as 



s{Ty 



1 - 2J=- J e-^^l^ sin2(2, r7)r/2d?7 , (4) 



where sin^ 6'(2;, r^) = with z = (3 / R{T) and e{r]) = r]/ cos2h{ri). 

We can obtain S(T) = R(T = 0)/R(T) if we know the temperature dependence of 
the order parameter on the left hand side of Eq.(^). As there are no phenomenological 
inputs for this, we shall consider the results from chiral perturbation theory (CHPT) which 
is expected to be valid at least for small temperatures. For higher temperatures near the 
critical temperature, lattice simulations seem to yield the universal behaviour with a large 
correlation length associated with a second order phase transition for two flavor massless 
QCD. We shall use such a critical behaviour to consider the temperature dependence of the 
order parameter near the critical temperature. 



We quote here the results of CHPT obtained by Gerber and Leutwyler ||T3[. The con- 
densate ratio at temperatures small compared to the pion mass is given as 



{(1Q)t _ 1 ^ c 



{qq)T=o F"^ 
where the functions h are defined as 

h = H\fi)/{27c') h[ = -H%f^)/{47c'T') 



(5) 



h = 3ho[hQ + fi'^hi] h = 3/;.q[/io + /i^/^-i] + 3/;,o[^o + ^i/^^ + /^^^i] 



(6) 



with = Mt^/T. Also 6e// 

3 / 1 35m^ 



0.6T 



^J-^l and a 



2a 



+ 



The constant c ^ 0.9 and F^/F = 1.057 ± 0.012 with F^ = 93MeV 



p!3| . The constants a and 6 are related to the S-wave and D-wave tt-tt scattering lengths 



respectively |]T3[. Finally the functions H^{^) are given as |]n 



VP+7^gAA2+^ 



(7) 

gV - ' — 1 

We have extracted the temperature dependence of the condensate as in Eq.(|) for low 
temperatures. For temperatures close to Tc, the critical behaviour is that of 0(4) spin 
model in three dimension and has also been seen in lattice QCD simulations [1^. The 
order parameter here is given as ^^^^^^^^ = (1 ~ ^)^! where (5 = 0.39 0. We have taken 
Tc = 150 MeV 0. The two regions are joined smoothly and the result is shown in Fig. 1(a). 
This result is fitted with Eq. (^ to determine S{T) = R{0)/R{T), which is plotted in 
Fig. 1(b). We shall use it to calculate the quark propagator and the hadronic correlation 
functions. 

III. QUARK PROPAGATION IN THERMAL VACUUM 



In the calculation of correlators, quark propagators enter in a direct manner and hence 
it is instructive to study aspects of the interacting propagator in some detail 0. 

The equal time interacting quark Feynman propagator in the condensate vacuum is given 



as Sq^p(^x 



where, 



?/'^(a;), ?/'^(0) y, which at finite temperature reduces to 
S{x, T) = -77-^ / e*''-^ cos20[sin2/i - 7 ■kcos2h]dk. 



2 (27r)3 J 
i 7 ■ X 



[/l(x)-/2(x)] + 



1 h{x) 



/•oo / sin hx \ 

Ii{x) = k [ coskx ; cos29dk, 

Jo \ kx I 



(8) 
(9) 



(10) 



hix) = / k sinfcxcos2^e-^ (^^'^ dk, (11) 
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with X = \x\, k = \k\. 

The free massive propagator, which can be derived from S{x,T) by the substitutions 
sin2/(A;) = m^/e and cos2 f{k) = k/e, is given as 



So{mq,x,T) 



TTig {mgKi{mgx) - 2/5 (x)) - i {mgK2{mgx) + 2Ie{u 



(27r)2s 

°° , o „ . , X f°° k"^ ( , sin/cz' 



(13) 



where /5(x) = / — sin(A;x) sin^ ^(i/c, I^ix) = [ — (cos /ex ; ]sm^ 6dk. 

Jo e Jo e \ kx J 

Kiijnqx) and K2{mgx) are the first and second order modified Bessel functions of the second 
kind respectively. 

In Fig. 2 we plot the two components Tr S{x, T) and Tr (7 ■ x)S{x, T) of the propagator 
for massless interacting quarks given by Eq. (^ at T = MeV, T = 100 MeV and T = 
135 MeV, corresponding to the chirality flip and non-flip components considered by Shuryak 
and Verbaarschot [0. The normalisation is discussed in our earlier work ||^ 

To compare with the constituent quark models with an effective constituent mass, we 
have also plotted the behaviour of free massive quark propagator with masses 100 MeV, 
200 MeV and 300MeV. In the chirality flip part, the propagator in the condensate medium 
starts from zero, consistent with zero quark mass at small distances, attains a maximum 
value of about 250 MeV at a distance of about 0.9 fm and then falls off gradually. Further 
the interacting propagator overshoots the massive propagators after about 0.6 fm. We also 
see that with increasing temperature, the chirality flip component has a lower peak and the 
position of the peak shifts towards higher distances indicating the decrease of the dynamical 
mass with temperature. 

In the chirality non flip part, the interacting propagator starts from 1, again consistent 
with the behaviour expected from asymptotic freedom. But at larger separation it falls rather 
fast indicative of an effective mass of the order of 150 MeV. These features are qualitatively 



similar to that of the quark propagator at zero temperature |]6|JT7[], though quantitatively 
there are differences. Also, the non-flip component falls faster with increase of temperature. 

Clearly therefore, similar to the situation at zero temperature, whereas a constituent 
quark description is adequate to describe the behaviour of the chirality nonfiip part of the 
propagator, it is not so for the chirality flip part. 



IV. MESON CORRELATION FUNCTIONS 



In our earlier work, we noted that phenomenology of correlation functions necessitated 
introduction of irreducible four point structure (or fluctuations of the condensate fields) in 



vacuum 



11]. In fact, the meson correlation functions were different from square of the two 



point function (propagator) and the difference could be expressed in terms of the four point 
function. The expression of the meson correlation function at zero temperature defined in 



our earlier work |11] can be extended to finite temperatures as. 



R{x,T)=Tr S{x,T)r'S{-x,T)r +Tr{\ 



Ux)T s(-f)r 



I). 



(14) 



Where J{x) = 'iljl^{x)raj3ip'l^{x), is a generic meson current with F being a 4 x 4 matrix 
(1,75,7^ or 7^75), X is a four vector; a and /3 are spinor indices; i and j are flavour indices. 
The field S(x) is the condensate fluctuation field introduced in Ref. |]1T[ to include four point 
irreducible structures in QCD vacuum.. 

Thus, at finite temperature the correlator (Eq. 0) is now the square of the interacting 
equal time thermal propagator plus the four point contribution at finite temperature. The 
thermal quark propagator was obtained in the earlier section. We keep the structure of the 



fields S(a;, T) the same as for zero temperature |11 



X 



:f) + s 



a/3 



X] 



lA (7*7^) a/3 (^ijk 0^(^) + /i2 ^»I3 0(^) 



(15) 

(16) 



where the first term corresponds to vector fluctuations and the second to scalar, /ii and /X2 in 
the above equations are dimensional parameters which give the strength of the fluctuations 



and (f)^{x) and (f){x) are vector and scalar fields such that, with \Q) as the ground state of 
QCD, we have 

(Q|(/.^(f)0^(O)|fi) = S'^gvix); {nimmi^} = 9s{x) (17) 

At finite temperature, the functions gv and gs will be temperature dependent. We do not 
know how to calculate it except for a general property that the effect of the four point struc- 
ture should decrease with temperature. We take here a simple ansatz for the temperature 
dependence of gv and gs, 

^5,y(x,T)= (^^Ml_^ gs,v{x,T = 0) (18) 

Similar to calculations at zero temperature, we shall consider the ratio of the physical 
correlation function to that of massless noninteracting quarks at zero temperature given as 

Ro{x)=Tr [So{x)r' So{-x)r 

. The normalised correlation functions thus defined as 

c(f.r)^^ (19) 

are plotted in Figure 3 for the pseudoscalar and vector channels. 

As expected (on physical grounds) the amplitude of the correlator decreases with increas- 
ing temperature. The peak of the vector correlator shifts towards the right after T = CQT^. 
We might remind ourselves that the position of the peak of the correlator is inversely pro- 
portional to the mass of the particle in the relevant channel 0]. 

The spatial hadronic correlators have been used to extract the hadronic screening masses 
and widths at finite temperature 0. To extract the hadronic properties at finite tempera- 
ture, we use a phenomenological parameterisation as is usually done in sum rule calculations 
I8| , p!9[] . We may note here however, that the phenomenological inputs are not available at 



finite temperature. The correlators are parameterised with the mass, decay width and the 
coupling of the particle to the vacuum, all three parameters being temperature dependent. 
We first express the correlator in terms of spectral density function. 
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Rph{x) 



Jo 



(20) 



Then we use the following phenomenological parameterisation for the spectral density func- 
tion mjii, 



p^(s) = 3Aj5(s - M^) + 1^ tanh 



AT 



9{s-So)+T^Sp6{s) 



(21) 



3s 
8^ 



tanh 



4T 



e{s 



(22) 



where A is the coupling of the bound state to the current, M is mass of the bound state and sq 
is the threshold for continuum contributions. The last term in Eq. (^) is the scattering term 
for soft thermal dissociations (mainly through pions), which exists only at finite temperature 
P^ . This term is given as 



1 rbT 
So = lim — / duj 
IpI^o 2n Jo 



oo 

2 / ^2 



X n{ 



IpIx — 

2T ' 



■ \V\X + Co", 

2T ' 



(23) 



The derivation of the above expression is slightly tricky and we have given it in the appendix. 
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Following Ref. ||T^ we take Sp 

The mass, threshold and coupling are then extracted such that the correlators as obtained 
from Eq. ^ agree with the normalised correlation functions as calculated by us (Fig. ^ ||11|| . 
This is done for each temperature. The results are plotted in Fig. ^ for the pseudoscalar 
channel and in Fig. ^ for the vector channel. The results are also shown in Table |. 



V. SUMMARY AND CONCLUSIONS 

As can be seen from Fig. 3, with increase in temperature, the correlation functions have 
a lower peak indicating lack of correlations with temperature. In the vector channel the 
mass of the p meson appears to decrease beyond 120 MeV. The threshold for the continuum 
also decreases around the same temperature. The behaviour with temperature of these 



quantities is qualitatively similar to that found by Hatsuda et al 119|. We have also plotted 
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the temperature dependence of the couphng of the boundstate to the current which decreases 
with temperature but rather slowly as compared to mass or the threshold for the continuum. 
The temperature dependence of these parameters can be used to calculate the lepton pair 
production rate from p in the context of ultra relativistic heavy ion collision experiments to 
estimate vector meson mass shift in the medium. 

In the pseudoscalar channel the mass remains almost constant till the critical temperature 



whereas the thershold and the coupling decrase with the temparature |]20|. We may note here 
that in the pseudoscalar channel, the contribution to the correlation function mostly comes 
from the fluctuating fields and the temperature behaviour as taken in Eq.(0) essentially does 
not shift the position of the peak whereas the magnitude of the correlator decreases. That 
is reflected in the above behaviour of the parameters in the pseudoscalar channel. We may 
note here that similar behaviour of pion mass becoming almost insensitive to temperature 
below the critical temperature was also observed in Ref. where correlation functions were 
calculated in a QCD motivated effective theory namely the Nambu- Jona Lasinio model. 

We would like to add here that the present analysis will be valid for temperatures be- 
low the critical temperature. Above the critical temperature there have been calculations 
essentilly using finite temperature perturbative QCD in random phase approximations pi 



However, in the region above Tc, nonperurbative features have been seen to exist from stud- 
ies in lattice QCD simulations [0. In view of this, one possibly has to do a hard thermal 
loop calculation where a partial resummation is done . 
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APPENDIX A: 



Here we shall derive the scattering term Sp. This may be calculated by considering the 
imaginary part of the longitudinal correlator for space like four momenta and can be written 



as 



ImUoo 

Pi{uJ,p)= (Al) 



which is explicitly written as 



X 5icu - El + E2)5^''\p- + k2){n2 - n^). (A3) 



Here Ei = \j kf + ml ; E2 = \l k^ + ml and rii = n{Ei) is the Bose distribution function for 
pions. 



In general the expectation of a vector current with respect to a pion state is given as 23 



{n{ki)\J^\n{k2)) = {ki + k2)^GM (A4) 

where, p = ki ~ k2 and Gt^{p) is the pion form factor with Gt^{Q) = 1 Substituting this in 
Eq. (^^) and integrating over k2 we obtain 



Sp{u;,p) = 2 X / —^{GAp)?S{^ -Ei + E- 2){2E, - ij)\n2 - ni) (A5) 



with k2 = p — ki. Next, since the delta function above contribute to space like (p^ < 0) 
region we write it as 

6{uj -Ei + E2) = 2E26{{uj - Elf - El)e{-p^) 

To simplify further, we may change the integration over three momentum ki to the integra- 
tion over energy Ei and the angle cosO^^. Performing the integration over angles restricts 

1 2 
the lower limit of the enegry integral Ei as -Eimm = 7:{^ + IpI'^), where, f = (1 !!v^i/2_ 

Thus we have 
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1 roo 

ptiuj,p) = -bf / Glip)i2E, - uj^in, - n,) (A6) 
with E2 — El — uj. Next, defining the varible x through Ei = ^{uj + leads to 



1 r°° 
8^ Jv 



(A7) 



We shall consider the longitudinal form factor Sp{uj,'p) in a frame which is at rest with 
respect to the medium which implies that p — 0. In this limit the constraint Q < uo < 
also forces uj to approach zero. However the above integral becomes increasingly large as 
p — > such that the integrated quantity of Sp{uj,'p) within the phase space for uj remains 
finite. Thus we first integrate over this region with p finite and then take the limit p ^ 0. 
Thus let 



/ = lim / duj^pUuj,p) = 



(A8) 



so that pf{uj,p) effectively becomes a delta function. Thus the spectral density reduces to 

S. 



.2\ '-'P 



lim = 5{lo )- 

p^->0 /TT 



(A9) 



We also note that there arises no ambiguity from the pion form factor as G,r(p = 0) = 1. 
Now the integral I can be written as 



1 r\'P\^ 
I = — lim / du'^ ""^ 
Svr \p\-*oJo 



dx x^ [n{{\p\x - uj)/2T) - n{{\p\x + cu)/2T)] 



(AlO) 



We change the integration variables by putting u; = |p|A and x = ^1+ \p\2{i_x'^) ■ Hence 
the spectral density function can be written as 

xy dy 



/ = lim / d\ \ , 



.\p\x-iu. .\p\x + uj, 
n( — — — n( 



(All) 



2T ' ' 2T 

In the limit of \p\ 0, we may Taylor expand the difference of the distribution functions in 
the square bracket of eq.( |All| ) and have 



n{ — — — ) — ni 



2T 



2T 



2x\pW1- X^^dn 



dX 
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Substituting back in eq. ( |All|) and performing an integration by parts for dX integration we 
have 

I = — C dxT dyni )y. (A12) 

In the hmit of vanishing pion mas we have / = 27rT^/9 so that Sp = 
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FIGURES 



(q q>T/(q q)T=o r(o)/r(t) 
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T (in MeV) T (in MeV) 

FIG. 1. Figure (a) shows quark condensate at finite temperature normalised to that at zero 
temperature obtained from CHPT and Lattice. Figure (b) shows R(0)/R(T) as determined from 
Fig (a). 
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FIG. 2. The two components of the thermal quark propagator, (a) Tr(S) and (b) Tr[(7 • x)S\ 
versus the distance x(in fm). The three lines, dot, short dash- long dash and solid corresponds to 
massless quark interacting propagator S(x,T) at temperatures of 0, 100 and 135 MeV respectively. 
The three lines, short dashed, dot-short dashed and long dashed correspond to a massive free 
propagator with a mass of 100, 200 and 300 MeV, respectively at T=135 MeV. 
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FIG. 3. The ratio of the meson correlation functions at finite temperature to the correlation 

functions for noninteracting massless quarks at zero temperature — — ) — distance x 

Rq{x, T = 0) 

(in fm). The solid, dashed, dotted and dot-dashed lines correspond to temperatures T =0 MeV, 



T =130 MeV, T =140 MeV and T =148 MeV respectively. 
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FIG. 4. The temperature dependence of mass, threshold (So) and couphng for the pseudoscalar 



channel. Tc =150 MeV. The vertical lines represent the errors obtained while fitting 
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Vector Channel 
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FIG. 5. The temperature dependence of mass, threshold (So) and couphng for the vector chan- 



nel. Tc =150 MeV. The vertical lines represent the errors obtained while fitting 
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TABLES 



TABLE I. Fitted parameters 



CHANNEL Temp.(MeV) M (GeV) A ^(GeV) 

Vector 0.780 ± 0.005 (0.420 ± 0.041 GeV)^ 2.070 ± 0.035 

100 0.771 ± 0.001 (0.411 ± 0.062 GeV)^ 1.897 ± 0.033 

130 0.774 ± 0.001 (0.402 ± 0.061 GeV)^ 1.672 ± 0.027 

134 0.762 ± 0.001 (0.391 ± 0.061 GeV)^ 1.586 ± 0.026 

138 0.741 ± 0.001 (0.375 ± 0.061 GeV)^ 1.468 ± 0.024 

140 0.726 ± 0.001 (0.363 ± 0.062 GeV)^ 1.393 ± 0.024 

142 0.706 ± 0.001 (0.347 ± 0.063 GeV)^ 1.303 ± 0.024 

144 0.679 ± 0.001 (0.325 ± 0.064 GeV)^ 1.192 ± 0.024 

Pseudoscalar 0.136 ± 0.00014 (0.473 ± 0.021 GeV)^ 2.110 ± 0.152 

100 0.136 ± 0.00014 (0.457 ± 0.021 GeV)^ 1.974 ± 0.118 

130 0.136 ± 0.00014 (0.402 ± 0.018 GeV)^ 1.574 ± 0.052 

134 0.136 ± 0.00014 (0.380 ± 0.017 GeV)^ 1.445 ± 0.038 

138 0.136 ± 0.00014 (0.351 ± 0.016 GeV)^ 1.284 ± 0.025 

140 0.136 ± 0.00014 (0.332 ± 0.015 GeV)^ 1.190 ± 0.020 

142 0.136 ± 0.00014 (0.310 ± 0.014 GeV)^ 1.084 ± 0.014 

144 0.136 ± 0.00014 (0.282 ± 0.013 GeV)^ 0.965 ± 0.011 

148 0.136 ± 0.00014 (0.201 ± 0.009 GeV)^ 0.970 ± 0.020 

150 0.136 ± 0.00014 (0.131 ± 0.006 GeV)^ 0.777 ± 0.020 
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